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ABSTRACT 


kThe  indeterminacy  of  the  reduction  of  an  arbitrary  matrix  to  Jordan 
normal  form  is  used  to  study  the  Lur’e  canonical  equations  for  a  linear 
control  system  with  a  single  nonlinear  actuator.  The  Lur'e  transformation 
is  shown  to  be  essentially  a  scaling  operation  applied  to  an  arbitrary 
transformation  which  reduces  the  system  to  Jordan  canonical  form.  In  the 
case  of  nonlinear  elementary  divisors,  the  proper  choice  of  canonical 
variables  simplifies  the  calculation  of  stability  conditions. 
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AN  ALGEBRAIC  PROBLEM  IN  CONTROL  THEORY 

* 

Arthur  Wouk 

1.  In  the  series  of  papers  in  the  early  forties  summarized  in  [1],  [2], 

A.  I.  Lur'e  showed  how  to  utilize  Liapunov's  second  or  direct  method  in  the 
investigation  of  the  stability  of  linear  automatic  control  systems  with  a  single 
nonlinear  actuator.  His  approach  consists  of 

1°  the  transformation  of  the  original  system  of  differential 
equations  via  the  so-called  Lur'e  transformation  into 
canonical  coordinates  in  which  the  construction  of  the 
Liapunov  function  is  direct,  and 

2  °  the  conversion  of  the  differential  problem  into  a  purely  algebraic 
problem. 

We  will  be  concerned  here  with  the  questions  of  the  existence  and  con¬ 
struction  of  the  Lur'e  transformation. 

A  system  with  a  single  nonlinear  actuator  is  described  by  a  system  of 
differential  equations  (cf.  e.g.  ,  Hahn  [5],  p.  42): 

y'  =  Ay  +  up 

u'  =  f(s),  s=(b,y)-hu,  h>0  , 
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where  '  =  d/dt,  A  is  a  real  constant  n  x  n  matrix,  b  and  p  #  0  are  real 
constant  n-vectors,  (,}  denotes  scalar  product,  and  f  is  a  real  continuous 
nonlinear  function  such  that  sf(s)  >  0  .  Here  the  vector  y  represents  the 
physical  state  of  the  system,  the  scalar  u  represents  the  amount  of  control 
to  be  utilized  and  the  vector  p  represents  the  distribution  of  the  control  u 
among  the  state  variables.  An  error  signal  s  is  constructed  from  the  deviation 
of  y  from  a  plane  normal  to  the  vector  b  ;  modified  by  a  multiple  of  the  current 
control  u,  it  signals  an  actuator  whose  response  is  characterized  by  the  non¬ 
linear  function  f(s)  .  In  many  physical  systems,  the  ability  to  control  state 
variables  is  limited  to  those  state  variables  which  represent  physical  velocities; 
only  their  derivatives  may  be  changed  (Newton's  Laws  of  Motion.)  Thus  a 
physically  significant  problem  is  one  in  which  A,  p,  f  are  assumed  as  given; 
the  physically  interesting  problem  of  stability  is  that  of  selecting  a  vector  b 
such  that  the  solution  y  =  0  of  (1)  is  asymptotically  stable  in  the  large,  (cf. 
e.  g. ,  [4]  p.  7-8)  for  all  functions  f  out  of  a  suitable  class  of  functions. 

(This  means  that  all  solutions  of  (1)  tend  to  zero  as  t  -*■  °o  . )  it  should  be 
noted,  however,  that  in  many  electronic  systems,  all  the  variables  in  the  mathe¬ 
matical  representation  may  be  controlled. 

We  note  first  (compare  Hahn,  [5],  p.  43}  that  the  transformation  to  the 
desired  canonical  coordinates  is  accomplished  by  a  Lur'e  transformation 

y  =  Bx  +  uq  , 


which  converts  (1)  into  the  form 


ERRATA 


p.  7,  bottom  line  comute 


p.  8,  line  12 


r  T?  r 

(X-X^  1,  (X-X2Z,...,  (X-X^  q, 


commute 

ri  r?  r 

(X-Xj)  ,(X-X2)  ,...,(X-X  )  q, 


p. 14, line  6 


2 1  — -  2 1 

r.  +  2r  +  ...  +  r  +1  3,3,0J 

1  K 


‘r  +2r2+-2r3+...  2^+1"  3,3,0 


p.  14,  line  8 


p.  26,  first  diagram 


1  1  0 


0  1  0  ; 


0  1  0 


1  1 


p.  28,  line  3 


cnaracteristic 


characteristic 


p.  29,  line  5 
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(2) 


x'  =  Dx  +  mf(s}  , 

I  s'  =  (c,x)  -  hf(s)  . 


T  T  T 

Here  D  is  the  Jordan  normal  form  of  A,  c  =  B  A  b  (  denotes  transpose)  and 
$ 

m  is  a  vector  with  m  =  1  or  0,  i  =  1,  . . .  ,  n  .  It  is  easy  to  see  that  we 

require  B  to  be  a  nonsingular  matrix  such  that  B  ^  AB  =  D  and  that  Aq  =  -  p  , 

-1 

B  q  =  -m,  or 


(3) 


ABm  =  p  . 


The  choice  of  iru  as  either  1  or  0  is  a  normalization  for  the  convenience 
of  analysis  and  exposition.  What  is  important  is  the  vanishing  or  non-vanishing 
of  rm;  should  im  vanish,  then  the  corresponding  canonical  variable  x.  is 
not  directly  controlled  by  the  actuator.  (It  may,  however-,  be  indirectly  controlled.) 
Lur'e  [1]  (see  also  [2],  pp.  38-47)  gives  certain  algorithms,  valid  for  the  case  of 
a  matrix  with  n  distinct  characteristic  values,  in  determinantal  form,  for  what 
is,  in  effect,  a  construction  of  B  and  q  .  In  these  determinantal  algorithms, 
the  geometric  meaning  of  the  manipulations  is  not  clear,  and  the  construction 
fails  for  the  case  of  multiple  characteristic  values.  We  will  exhibit  a  procedure 
which  determines  the  possibility  of  construction  of  the  Lur'e  transformation  as  it 
constructs  it;  it  will  be  seen  that  the  algorithm  of  Lur'  e  valid  for  matrices  with  n  distinct 
* 

Hahn([5],  p.  43,  p.  138)  first  observed  that  it  is  not  always  possible  to  have 
mj  =  1,  i  =  1,  . . . ,  n;  i.  e. ,  some  zero  values  may  be  necessary  if  p  is  orthogonal 
to'  characteristic  vectors  of  A  . 


-4- 


#392 


characteristic  values  results  in  no  more  than  the  special  scaling  of  an  arbitrary 
B  which  transforms  A  to  Jordan  normal  form  so  as  to  set  rn.  =  1,  or  0  .  We 
will  show  that  the  question  of  which  canonical  variables  can  be  controlled 
(m.  =#  0}  can  be  resolved  by  knowing  any  matrix  B  which  reduces  A  to  Jordan 
normal  form. 

Recently,  Yakubovich  has  shown  that  the  transformation  to  canonical 

coordinates  is  unnecessary  from  the  point  of  view  of  theory  (cf.  [3]  pp.  75-103, 

in  particular  pp.  75-91  or  [8]  pp.  129-123  for  expositions  in  English  of  his  work.) 

That  is  to  say,  the  conditions  of  Lur'e  for  the  stability  of  the  controlled  system 

are  derived  by  matrix  theoretic  methods  which  do  not  directly  utilize  the  Jordan 

normal  form.  Hidden  in  the  theory  is  the  necessity  for  solving  the  matrix 
T 

equation  A  L  +  LA  =  -  C  for  any  prescribed  positive  definite  symmetric  matrix 
C  .  In  fact  this  is  a  very  difficult  matter  if  A  is  not  in  Jordan  normal  form. 

Thus,  at  the  point  of  application  it  remains  necessary  ,  (  cf.  e.  g.  LaSalle  and 
Lefschetz  [3]  pp.  91-103  for  the  case  where  this  form  is  diagonal), to  utilize  the 
transformation  to  canonical  coordinates  in  order  to  establish  computable  constraints 
on  the  control  variables  p,  b  and  h  of  (1)  which  lead  to  desired  stability 
properties.  (In  this  regard,  see  also  [4].)  Thus  it  remains  necessary  to  determine 
which  sets  of  canonical  variables  may  be  achieved  by  this  reduction,  what 
effects,  if  any,  this  indeterminacy  in  the  reduction  has  on  stability,  and  which 
canonical  variables  permit  easiest  calculation  of  stability  criteria.. 

Lastly  in  an  appendix  we  will  relate  these  results  to  those  of  Kriuchkov  [9] 
which  construct  a  "canonical  form",  other  than  the  Lur'e  form,  which  need  not 
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preserve  the  stability  properties  of  the  original  system.  We  will  see  that  the 
Lur'e  form  is,  in  suitable  modification,  always  achievable,  and  that  stability 
is  always  preserved. 

2.  We  note  first  that  if  we  set  p  =  Bp,  then  (3)  becomes  Dm  =  p  .  Thus, 
for  a  given  B,  the  decomposition  of  p  along  the  basis  consisting  of  the  column 
vectors  of  B,  b^,  ...,  b^,  determines  the  possibility  of  finding  m  .  If  m 
exists,  it  is  unique  up  to  a  solution  of  the  homogeneous  equation  Dm  =  0,  i.  e. , 
to  an  arbitrary  vector  in  the  null  space  of  D  while  p  must  be  orthogonal  to 
the  null  space  of  D  .  Thus  p  must  be  orthogonal  to  the  null  space  of  A  . 
Henceforth  this  is  always  assumed. 

Next  we  note  that  if  \  #  0  is  a  characteristic  value  of  A(and  D)  possessing 
two  or  more  linearly  independent  characteristics  vectors,  then  there  exist  two 
columns  of  B,  b,  .  b.  satisfying  Ab.  =  X.  b.  ,  j  =  1,  2  .  If  p  and  p, 

7  i  7  i  i  i  7  7  ri  are 

1  2  ]  ]•  X1  *2 

both  nonzero,  than  we  must  have  p,  =  p.  .  To  see  this,,  note  that  Dm  =  p 

h  12 

implies  X  ,m,  =  p,  .  Since  X.,  p,  and  p.  are  not  zero,  then  we  must  have 

11.  1 1  1*  li  1« 

]  J  12 

m.  =  1,  j  =  1,  2  and  p,  =  p.  =  X  .  .  This,  however,  is  only  an  apparent  restriction, 

X1  *2  1 

for  we  can  always  rescale  p  by  renormalizing  the  characteristic  vectors  b.  , 

X1 

b,  ,  ...,  so  that  p,  =  (b.  ,  p)  is  a  constant  for  all  characteristic  vectors 
x2  xj  Xj 

belonging  to  a  given  nonzero  characteristic  value.  Note  that  this  entails,  with 

therenormalization  of  b.  ,  the  renormalization  of  the  other  columns  of  B  which 

j 

are  in  the  same  cyclic  invariant  subspace  as  b, 

j 

On  a  given  cyclic  invariant  subspace  corresponding  to  an  elementary  divisor 

k. 

(k  ~\  .)  ,  D  has  the  representation 
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(4) 


V  . 

l 


1 

X .  1 


V  . 
1 


1 

X 


If  the  relevant  components  of  m  and  p  are  indexed  by  j  ,  . ..,  j,,  we  find 
that  on  this  subspace  {3}  is  equivalent  to 


X.  ,m.  +  m.  =  p. 

1  h  h  J1 


\  ,  m ,  +  m ,  =  p, 

1  J2  J3  h 


X.  ,m.  +  m,  =  p, 

1  jk-l  Jk  Jk-1 


X  ,m. 

1  Jk 


From  these  equations,  it  becomes  clear  that  (3)  is  not  in  general  solvable  for 

im  =  1  or  0,  and  the  construction  of  the  transformation  is  not  immediate.  Further 

it  seems  that  renormalization  of  the  b,  might  not  effect  equality  of  other  p  's 

Jk  1 

belonging  to  distinct  invariant  subspaces  corresponding  to  the  same  X  .  It 

i 

might  appear  from  this  that  the  restriction  of  m,  to  the  values  1  or  0  is  at  fault 
here.  If  im  were  allowed  arbitrary  values,  then  the  numerical  value  of  the  p, 
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would  enter  into  the  determination  of  the  m.  ;  we  shall  see  below  that  this  is 

1 

misleading:  only  the  vanishing  and  non-vanishing  of  some  of  the  p  are 

l 

significant  to  this  determination. 

3.  The  starting  point  of  our  procedure  is  the  observation  that  the 

transformation  to  Jordan  form  of  A  is  unique  only  to  an  arbitrary  nonsingular 

-.1  -1 

matrix  V  which  commutes  with  A,  for  then  (VB)  -A(VB)  =B  AB  =  D .  Hence- 

forth  we  assume  that  we  know  a  matrix  B  which  reduces  A  to  D.  The  construction 

of  such  matrices  B  is  well  known  ( cf.  e.g.,  [6]  p.  66-69  or  [ 7 ] p.  159  ff .  )  .  We  can 

restate  our  problem  as  the  determination  of  a  matrix  V  which  commutes  with  A  such 

that  under  the  transformation  y  =  VBz  +  uq  the  system  (  1)  goes  over  into  (  2)  where 
T  T  T 

c  =  B  V  A  b,  and  in  addition 

(  5 )  AVBm  =  p . 

Hence  q  =-VBm.  (We  omit  components  in  the  null  space  of  A  . ) 

In  terms  of  the  given  nonsingular  matrix  B,  p  can  be  decomposed 
uniquely  as 


P  =  E  =  BP 
i=i 

for  a  unique  vector  p,  with  components  ,  ...  }  p  .  Among  the  matrices  which 
comute  with  A  are  the  polynomials  in  A,  say 
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n  -  1 

g(a)  =  ^  ^A1  • 

i=  1 

Unfortunately,  if  any  two  elementary  divisors  of  A  have  a  nonlinear  common 
factor,  a  suitable  V  may  be  found  in  this  way  only  for  restricted  classes  of 
vectors  p  ( see  below) .  Instead,  we  observe  (cf.  e«g.  [6],  pp.  143  ff.  or  [7] 
pp.  220-223)  that  the  set  of  matrices  V  which  commute  with  A  is  generated  by 
V  =  B  UB  1  where  U  runs  over  the  set  of  matrices  which  commute  with  D.  Then 

(5)  can  be  rewritten  as  B  ^AVBm  =  p  or 

(6)  DUm  =  UDm  =  p. 


4.  In  order  to  construct  the  desired  nonsingular  matrix  V,  we  must 
study  the  block  structure  of  U,  which  is  dependent  upon  the  elementary  divisors 
of  D .  If  these  are 

r]  r?  rrr 

(X  -xx)  \  (X  ... ,  (X  -xq)  q, 

then  ([6]p.l44or  [7]p. 218)  U  is  decomposed  into  blocks  TJ  ,  or(3  =  1 ,  ...  ,  q 

“P 

where  U  „  is  a  rectangular  matrix  with  r  rows  and  r_  columns.  Further, 

up  a  p  3 

U  _  is  a  zero  matrix  if  X  +  X.,  while  it  is  "upper  triangular"  if  X  =  . 

ap  up  a  p 

( "Upper  triangular"  for  an  m  row,  n  column  matrix  M  means  M.,  =  0  if 
j  <  i  +  max  ( 0 ,  n  -  m) ;  for  square  matrices  this  is  the  usual  definition. )  Lastly, 
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the  elements  along  each  diagonal  of  U  .  are  constant.  Thus  each  block  contains 

ap 

min  (r  ,  r  )  arbitrary  parameters.  It  would  therefore  appear  that  we  have  sufficient 

a  p 

freedom  to  restrict  the  class  of  matrices  U  among  which  we  look  for  our  solution, 

possibly  to  upper  triangular  matrices  ( =  the  zero  matrix  if  a  >  (3. )  It  will  be 

seen  below,  however,  that  this  can  be  done  only  in  special  circumstances;  this 

will  explain  in  part  our  inability  to  find  U  among  the  matrix  polynomials  in  D . 

We  will  prove  first  two  lemmas,  one  guaranteeing  the  solvability  of 

(6)  for  arbitrary  p,  the  other  guaranteeing  the  nonsingularity  of  U. 

It  is  evident  from  the  description  that  we  can  restrict  our  attention  to  the 

invariant  subspaces  associated  with  a  given  nonzero  characteristic  value,  say  0 

X^,  since  U  and  DU  are  reduced  to  block  diagonal  form  with  respect  to  these 

subspaces.  Let  W  denote  the  block  in  U  which  corresponds  to  X^.  We  rewrite 

rl  rk 

the  form  of  (6)  on  such  a  subspace  as  follows.  Let  (X  -X^)  ,  ...  ,  (X  -X^)  , 

r^  _>  r  >  .  .  .  >  r^  >  1  be  the  elementary  divisors  associated  with  X^  .  Let  the 

index  pair  [i,  t]  denote  the  index. 

i  -  1 

Z  rf3+t 

(3=1 

where  t  =  1 ,  ...  ,  r. ,  i  =  1 ,  2 ,  ...  ,  k .  ( The  sum  is  taken  as  zero  when  i  =  1  . ) 

Let  (U  a)  .  .  =  z  .  .  .  .  ,  _  .  , 

ap  l,  j  a,  (3,  i  -  j ,  i  =  1 ,  ...  r  ,  j  =  1 ,  ...  r  ,  a,  p  =  1 ,  ...  ,  k  .  Since 

O'  p 

U«(3  is  upper  triangular,  it  follows  that  z  a  .  .  must  vanish  if  j  <i+  max(0,  r„  - r  ). 

a,  p,  i  -  j  p  a 

On  the  other  hand,  the  matrix  product  DU  =  UD  has  the  same  block  structure  as 
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U,  with  elements  z1  „  ,  ,  =  X  z  Q  .  .  +  z  (We  extend  the 

5  a,  (3,  l-j  1  a,  p,i-j  P,i-J+1 

definition  of  z  n  .  ,  ,  ......  ,  ^  .  ,  ,  .  ,  ,  , 

a,  p}  i-j  by  setting  it  equal  to  zero  if  l  or  j  are  outside  their 

respective  ranges  1,  ...  ,  r  and  1.  ...  r„  )  .  Let  W*  be  the  block  in  DU  which 

a  p 

corresponds  to  W.  Then  (6)  becomes 


(7) 


k 

E 

p=i 


E 


“  p  r 


i  =  1 


>  Va=1> 


If  we  take 

r.  +  ...  +  r. 

1  k 


m,fl  as  given,  and  hopefully  =  i, 

L  Pj  J  J 


equations  in  the  unknowns  z1 


a,  P,  i-j 


then  (7)  is  to  be  viewed  as 
There  are 


k  k 

m-Tj  Tj  min(r  ,  r)  =  r  +3r  +...+( 2k- 1)  r 
as  1  p  =  l  P  1  ^  1 


such  unknowns.  Thus  we  have  more  unknowns  than  equations.  Solutions  exist 
if  the  equations  are  consistent.  The  consistency  conditions  may  be  expressed 
in  terms  of  an  auxiliary  matrix  M  which  we  determine  thus:  enumerate  the  z1 
variables  so  that 


z 


t 

1 


-  r 


1 


+  1 


2,  0’ 


+  r 


1 


r 


2 
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;1'r  + ...  +rk_L+l  =  2,1 1}  k.  O’  -  ’  Z,ri+...  +rk“  Z\,k,l-vk  ’ 


Z*r^  + ...  +^k  +  l  “  Z  2,  lj  l-(r^  - r^) *  ’  Z  1^  +  ...  +rk  +  r2_  2  2#1>1“r1 


ZV  +  ...  +  r  +  r  +  1  =  *2,2,0’  * 

1  k  2 


That  is,  we  start  with  the  upper  left  corner  of  the  W*  matrix  and 

enumerate  from  left  to  right  in  each  row,  those  elements  zl  ,  in  the  l*st 

r> 

r  +1*  st,  ...  ,  r,  +  ...  +r,  ,  +  Pst  rows  for  which  j 1  >  1  +  max  (0,  r„  -  r  ) . 

1  1  k-1  —  pa 

Then  Mz1  »  p  where  M  is  a  rectangular  matrix  of  the  following  block  form: 


where  is  a  matrix  of  r  rows,  ir^  +  r^  +  ...  +  rk  columns.  If  we 
subdivide  Mk  into  blocks  of  r^,  ...  ,  r^,  r^j,  •••  >  rk  columns,  then  the  j‘th 
block  has  the  form  of  the  following  r,  x  r.  matrix 
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i 

i 

i 


if  j  <  i,  while  it  has  the  form 


if  j  >  i. 
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From  this  structure,  it  follows  that  the  assignment  mrR  ,  =1  provides 

LP>  tJ 

a  matrix  M  which  is  of  maximal  rank,  since  its  rows  constitute  an  independent 
set  of  r  +  ...  +  r  vectors  in  an  m-dimensional  space.  In  fact,  it  is  sufficient 

J-  K 

for  this  if  mr  •,  4-  0;  then  the  l‘st  block  of  M.  is  a  square  nonsingular 
’  rl  J 

matrix.  On  the  other  hand,  mr,  ,=  0  makes  the  rank  of  M  less  than  r, +  ... +r 

[i,  fji  ik 

so  that  the  existence  of  DU  solving  (6)  is  not  guaranteed  for  arbitrary  vectors  p. 

We  have  proven  the  following: 


Lemma  I:  mr,  ,  =  1 ,  is  necessary  and  sufficient  for  the  solvability  of  (7) 

_[  !  >  rj[]  - 

for  arbitrary  p . 

k 

There  are  Y  (2l-2)r.  independent  parameters  in  the  W1  matrix  whose  assignment 
j  =  l  J 

.leads  to  a  solution  of  (7)  with  mr  .  =  1 ,  p  =  1 ,  ...  ,  k,  j  =  1 ,  ...  ,  rR .  If,  as 
would  be  desirable  for  computational  purposes,  the  matrix  U  is  sought  among  the 
class  of  upper  triangular  matrices,  the  following  changes  must  be  made  in  the 
preceding  analysis.  W  and  W1  are  now  upper  triangular  matrices  and  ( 7)  specializes 
to 


( 71) 


z 


Pj  i-jm[P>  J] 


=  Pr 


P -« j- 1 


i  a  1 


>r  ,  a  -  1 , 

a 


a,  p,  i-j 


The  number  of  unknowns  z1 


is  now 
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m  =  f  l  min  (r  ,  r  )  =  r  +  2r  +  ...  +  kr 
a  p  I  £  k 


ors  1  (3  =  1 


which  still  exceeds  the  number  of  equations.  We  proceed  as  before  in  the  enumeration 
of  the  elements  of  zl ,  i.  e.  from  left  to  right  in  W* ,  but  now 


11,1,0’ 


+  r^  +  ...  +r^+  1  2,  2,  0* 


r_  +  2r  +  . , .  +  r.  +  1  3,3^0‘ 

12  k 


M.  is  now  a  rectangular  matrix  of  r.  rows,  r.  +  r.,  ,  +  ...  +  r,  columns  whose  j:th 
i  l  l  l+l  k 

column  block  is  given  only  by  the  second  of  the  forms  above,  that  for  j  =  i+  p  -1>  i, 

The  assignment  m^  t]  ~  1  Provic^es  a  matrix  M  of  maximal  rank,  as  before. 

However,  in  this  case  mr .  ,  *  1  makes  the  first  column  block  of  M.  a  non- 

[Ja  r  ]  j 

singular  square  matrix,  while  if  any  m^  f  ^  =  0,  the  corresponding  M,  has  a 
row  of  zeros.  We  have  then 


Lemma  I1:  _m^  r  j  =  1 ,  j  =  1 ,  ...  ,  k  is  necessary  and  sufficient  for  the  solvability 
of  { 71}  for  arbitrary  p . 
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q. 

There  are  Yi  (j  -l)r.  independent  parameters  in  the  W*  matrix  whose 
j  =  l  J 

assignment  leads  to  a  solution  of  (71)  with  mrR  =  1 ,  j=  1,  ...  ,  r  ,  p  =  l,  ...  ,k. 

LP>  Jj  P 


5.  We  wish  next  to  guarantee  the  nonsingularity  of  U .  Again,  owing 

to  the  block  structure  of  U  and  DU,  we  work  on  the  invariant  subspace  associated 

with  X.  4-  0;  there  the  nonsingularity  of  DU  and  U  are  equivalent.  The 

singularity  of  DU  or  W1  there  is  equivalent  to  the  existence  of  a  set  of  numbers 

6  r  o  ,, ,  (3  =  1 ,  ... ,  k ,  j  =  1 ,  ...  ,  rft ,  not  all  zero,  such  that 
l  P  s  )  J  ■  P 


(8) 


k  ^ 

pL1  j”!  Z  P>  i“j  6[P>  l]  "  °J 


Set  i  =  r^,  a=  1,  ...  ,  k,  and  consider  the  resulting  k  equations 


(9) 


Z 


6[P.rp]’0’ 


or  =  1  ,  ,  k 


If  the  determinant  of  the  auxiliary  k  x  k  matrix  (  zl  „  )  is  not  zero,  then 

a.  p,  r  -r„ 

'  a  p 


we  must  have  5 


[ P ,  rp]  =  0,  P  -  1,  ...  ,k,  and  ( 8)  can  be  replaced  by 


ra-l 


(10) 


Z  Z 

P=1  )=1 


P,  i-j  [P>  j] 


6rn  n  “  °»  i=  1  ,  ...  r  -1,  a- 1 ,  ...  k  . 


Now  set  i  =  r  -1,  a  =  1 ,  ...  ,  k.  There  results  the  system  of  equations 
a 
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js. 

^  2,1  a,  6.  r  -r  6  [3,  r  -  1]  ~  °* 

p  - 1  5  p>  a  p  1,  5  rp  J 


os  lj  ...  j  k-1  . 


The  nonvanishing  of  the  determinant  of  the  auxiliary  matrix  again  implies 

6  n  =  0,  ps  1,  ...  ,  k.  We  may  proceed  in  this  manner  for  exactly  r 

L p >  rp . "  1 J  k 

steps.  Suppose  that  r,  >  r  if  is  1,  s  while  r.  =  r  is  is  s  +  1,  ...  ,  k. 

X  K  X  K 

Then  it  follows  that  z*  .  =  0  if  a  <  s  and  P  >  s  +  1  .  Then 

«•  r«  -  rp 


det  (z1 

1  <  a,  p  <  k  a, 


)  =  det,  ,  A  (z1  _ 

■rp  l<^P<sv  «»P»ra-rp 


)  x  detg  +  ^  <  a  p  <  a  p  r  _  r  )  • 
—  —  a  p 


Then;  at  the  end  of  r^  steps,  the  successors  to  equations  ( 10)  become 


(10*) 


s  Vrk 

Z  Z 

p=i  j=i 


ZjPji-j  [Pjj] 


n  = 


i  =  l,  ...  ,  ra  -rk,or  =  l,...  ,s. 


Now  only  the  s  x  s*th  principle  minor  of  the  auxiliary  matrix  remains.  This  is 
non-singular  if  and  only  if  the  auxiliary  matrix  is  non-singular.  We  may  proceed 
in  this  manner  to  show  that  if  det  ( z*  „  )  #  0  ,  then.  W5  is  nonsingular. 

a)  Pi  -  rp 

On  the  other  hand,  if  det  (z1  R  )  =  0  then  (9)  has  a  nonzero  solution; 

a  p 

then  we  can,  by  elementary  column  operations  on  W* ,  construct  a  matrix  k  of 
whose  rows  vanish,  so  that  the  matrix  is  singular.  Since  elementary  column  operations 
are  themselves  nonsingular,  this  proves 
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Lemma  II: 


(11) 


det  ( zJ 


P>  r  -r. 


)  *  0 


is  necessary  and  sufficient  for  the  nonsingularity  of  W1 


If  all  r  are  distinct,  so  that  ( z1  a  ) 

matrix,  or  else  if  we  seek  an  upper  triangular  matrix,  U  , 


is  a  lower  triangular 
condition  ( 1 1 )  becomes 


(lli) 


k 

n 

a=l 


z1 

a, 


a,  0 


#  0  . 


6.  Next  we  establish  the  conditions  for  the  existence  and  nonsingularity 

of  the  matrix  U .  We  may  still  restrict  attention  to  the  block  W1  in  DU 

corresponding  to  a  single  nonzero  characteristic  value  .  In  (7),  set 

i  =  r  ,  a  =  1 ,  ,  k  .  This  yields  the  set  of  equations 

a 

k 


(12) 


Z 

(3  =  1 


a,  0,  r  -  r  a 
7  a  p 


m 


[P,rp] 


=  P 


a  =  1  ,  ...  ,  k  . 


If  p,  ,  =  o ,  a  =  1 ,  ...  ,  k,  then  it  is  impossible  to  have  any 
L  > 

mr  n  +  0  and  simutaneously,  a  nonsingular  z*  -  matrix  (or,  by  Lemma  II  , 
[a,  r  J 

L 


a  nonsingular  matrix  W1 . )  Thus  we  must  set  mr  =  0 ,  a  =  1 ,  ... ,  k  and 

I  a,  r  1 
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can  choose  z1  „  at  will  so  as  to  satisfy  ( 11)  . 

Next  consider  the  situation  where  p;  ,  ?=  0.  Suppose  first  that 

[a,  r  ] 

r  >  r„  .  Then  zJ,  =  0,  a,  =  2,  ...  ,  k;  if  pr.  ■, # 0 ,  then  (  12)  may  be 

1  2  1,  o,r1  -rk 

solved  for  a  nonsingular  lower  triangular  z1  -  matrix  and  a  vector  mrQ  n  -  1 , 

[ P>  rp] 


P  a  1,  ....  k.  We  need  only  set  z1 ,  ,  .  =  pr ,  n 

1  j  1>  0  [If  rj] 


and 


1  1  +  Z*  r,  -  Pr  1) 

^  1,ra"ri  01,01,0 


^  — *  1  j  «  • «  j  k  j 


while  all  other  z*  ss  0  .  This  can  always  be  done  so  that  both  (11)  and 

P’  ra“rp 

( ll1)  are  satisfied  and  the  matrix  W*  is  nonsingular. 

If  r  s  r_  «  ...  ar  j  s  <  k,  and  at  least  one  of  pr  n  ^0, 1  <  a  <  r, 
12  s'  -  [«»ra]  -  - 

then  we  may  reorder  the  basis  so  that  p^  ^  ^  is  again  not  zero,  and  proceed 

as  above.  It  is  also  possible  to  have  mr,  n  -  1,  mrra.  «  0  or  1  at  will, 
- e -  [l,rx]  [P,rp]  “  - 

(3  ts  2,  ...  ,  k,  and  preserve  the  non-singularity  of  WJ .  Thus  to  have  mrra  ^C), 

LPj  r-pj 

set  z1  „  ,  -  Pro  u  z*  o  #  0 . 

P>  1,  rp-rx  F[P,  rp]  P,  P,  0 


If,  however,  Pr  ,  a  0,  as  1,  ...  ,  s,  then  we  must  set  mr  -  a  o, 
«,r  [ a,  r  ]  ’ 

«  ss  1,  ...  ,  s  ..  Otherwise,  the  matrix  whose  elements  are  z1  R  a,  p,  =  1,  ...  , 

“»  p>  Vr(3’ 

would  be  singular.  Since  the  other  elements  of  the  first  s  rows  of  the  zJ  matrix 
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are  zero,  this  implies  the  singularity  of  the  z*  matrix,  and  by  Lemma  II,  the 
singularity  of  the  W1  matrix.  When  mr  =  pr  ,  s  0,  as  1,  ,  s, 

la>rJ 

the  quantities  z1  p  a,  p,  a  1 ,  ...  ,  s  may  be  chosen  as  an  arbitrary, 

a’p»ra-rp* 

nonsingular  s  x  s  matrix.  The  question  of  the  non  singularity  of  the  z*  matrix 
is  now  reduced  to  that  of  the  nonsingularity  of  the  k-s  x  k-s  matrix  z1  _ 

^P*ra-rp 

with  a,  p  a  r  +  1 ,  ...  ,  k.  The  argument  above  applies  without  exception  to  this 
submatrix. 

Now  suppose  that  m,fl  ,,  p  a  1 ,  ...  ,  k,  and  the  nonsingular 

LP>rpJ 

auxiliary  k  x  k  matrix  ( z1  )  have  been  chosen.  The  remaining 

0!m  r  •*  IT  n 

7  a  p 

equations  of  (7)  can  be  written  as 


(13) 


kJ  rp-l 

i?i  j?i  i“l>~.,ra-l,...,k, 


where  k'  _<  k  is  such  that  r  a  1 ,  k'  <  a  <  k .  With  appropriate  identification  of 
terms,  this  system  is  of  the  same  form  as  (7).  Hence  Lemma  I  may  be  applied  to 
(13).  We  conclude  that  the  assignment 


m 


[P^p-1] 


a  1 


P« 


k' 


is  necessary  and  sufficient  for  the  solvability  of  (13)  for  any  vector.  The 
remaining  elements  of  the  m -vector  may  be  chosen  as  0  or  1  at  will. 


-20- 


#392 


We  summarize  these  results  as  follows: 

Let  X  #  0  be  a  characteristic  value  of  A  with  elementary  divisors 

r  r 

1  k 

(X-X  )  ,  ...  ,  (X-X  )  ,  r  >  r„>  ...  >  r  >  1  and  characteristic  vectors 

1  1  -L  l  K 

b,  ,  ...  ,  b,  .  Let  ( b,  ,  p)  =  0  ,  i=  1 ,  ...  ,  s  -  1  while  ( b_  ,  p)  +  0  .  Let 

1  K  1  S  " 

r  .  >  r  =  r  =  r  =  ...  =  r  >  r  ,  ,  .  Then  the  assignment 

or-1  a  a  +  1  s  t  t+1  - 

s  s  s  s  s 


mr  ,  =  0, 

a,r 
L  or 


a 


a  -  1 
s 


m  r  1  =  1} 

[a  , r  ]  5 

s  a 

s 


assures  the  existence  and  permits  the  calculation  of  a  nonsingular. matrix  U 
which  commutes  with  D  and  satisfies  (6)  .  The  remaining  elements  of  m  may  be 
taken  as  1  or  0  at  will. 


The  choice  of  m 


fa  ,  r  1  in  the  above  statement  is  not  essential:  any 
L  s  «  J 
s 


one  of  mr  t  could  be  chosen  to  be  1 ,  a  <  a  <  t  and  the  statement  remains 

a.  r  I  s  —  —  s 

L  aJ 

true.  Further,  at  least  one  of  them  must  be  so  chosen. 

The  matrix  U  can  be  chosen  among  the  upper  triangular  matrices  if  and 
only  if  ( b^,  p)  #  0  for,  at  least  one  characteristic  vector  b^  belonging  to  each  set 
of  elementary  divisors  of  given  exponent.  To  see  this,  consider  again  equations 


(71)-  Here  z*  „  =  0  if  r  ^  r„ .  Repeating  the  preceding  analysis  leads 

or,  (3,  r2-rp  a  p 


to  this  form  of  the  assertion,  using  Lemma  I*  in  the  place  of  Lemma  I.  Lastly, 
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if  the  elementary  divisors  have  no  nonlinear  common  factors,  then  U  may  be 
realized  as  a  polynomial  g(D)  ,  after  renormalization  of  the  corresponding 
characteristic  vectors.  To  see  this,  note  that  this  implies  that  U  may  be 
taken  as  a  block  diagonal  matrix,  each  block  of  which  is  precisely  of  the  form 
which  is  realizable  by  a  suitable  interpolating  polynomial  in  D  ( see,  for  instance, 
[6],pp.  62-63  or  [7],  p.  100  .)  If  the  elementary  divisors  are  all  linear,  then 

D  is  diagonal  and  U  is  determined  by  (6)  as  a  diagonal  matrix,  used  only 
for  scaling  purposes. 

Thus  the  particular  transformations  described  by  Lurle  are  seen  to 
effect  only  the  scaling  of  rm,  when  they  are  correct,  i.  e. ,  in  the  case  that  A 
has  n  distinct  characteristic  roots. 

The  preceding  treatment  leads  to  complex  transformations,  if  the 

matrix  A  has  complex  characteristic  values.  Since  Ax  =  Xx+y  if  and  only  if 

Ax'-  X  ‘x’  +  y"  ( *  denotes  complex  conjugate)  when  A  is  real,  the  equivalent 

real  transformations  can  be  constructed  which  lead  to  the  real  Jordan  normal  form, 

in  the  usual  manner.  The  reality  of  p  and  b  implies  that  all  operations  associated 

$ 

with  the  characteristic  value  X  are  exactly  parallelled  in  the  treatment  of  X  , 
and  the  final  form  (2)  can  be  chosen  with  the  control  f(  s)  present  ( m.  =  1)  in 
precisely  one  pair  of  equations  belonging  to  each  pair  of  characteristic  values 
X,X\  for  which  (p,  b  )  ^0  for  some  non-zero  characteristic  vector  b^ 


belonging  to  X . 
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7.  The  indeterminacy  in  the  form  of  m  allows  us  to  raise  two  questions. 

Firstly,  does  this  indeterminacy  lead  to  different  stability  criteria?  Secondly, 
does  this  indeterminacy  permit  simplified  calculation  of  conditions  which 
guarantee  absolute  stability?  The  first  question  is  answered  as  follows.  Suppose 
that  a  matrix  B  has  been  calculated  and  the  equations  are  in  a  real  Jordan  normal 
form  (  2) ;  let  another  canonical  form 


(2') 


x1  =  Dx  +  mf(s) 

< 

,  S*  =  (C,  x)  -  hf  (s) 


be  achieved  by  the  matrix  U  where  U  commutes  with  D  and  is  non-singular. 

—  -1  —  t 

Then  m  =  U  m  and  c  =  U  c.  The  Lur:e  special  Liapunov  function  implies 
absolute  stability  if  and  only  if 


(14) 


h  >  ( A_1[Lm  +  £c] ,  [Lm  +  |c] ) 


where  A  is  some  real  symmetric  positive  definite  n  x  n  matrix  and  L  is  the 
unique  real  symmetric  positive  definite  solution  of  the  matrix,  equation  D^L  +LD=  -  A 
(cf  [3],  p.  85  or  [8],  p.  122.)  Then 
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while 

DT(  UTLU)  +  (  UTLU)  D  =  -  utau  . 

From  this  we  conclude  that  a  triple  (h,  m,  c)  makes  (2)  absolutely  stable  if  and 
only  if  the  triple  (h,  m,  c)  makes  (2l)  absolutely  stable.  This  answers  the  first 
question  negatively.  The  second  question  can  be  answered  affirmatively.  Notice 
first  that  the  stability  criterion  (14)  depends  upon  the  calculation  of  the  solution  of 

(15)  DTL  +  LD  =  -  A  ; 


this  is  not  readily  calculable  for  arbitrary  positive  definite  matrices.  Even  for 
3x3  matrices  the  computations  are  too  cumbersome  unless,  as  is  customary,  we 
make  A  a  diagonal  matrix.  This  yields  sufficient  conditions  for  (14)  which  in 
turn  is  only  sufficient  for  absolute  stability.  If  we  solve  ( 14)  for  all  sets  of  h 
and  c  which  satisfy  ( 14)  for  a  given  m  and  some  positive  definite  diagonal 
matrix  A ,  we  obtain  stability  criteria  which  depend  upon  the  form  of  m .  If  X  is 
a  characteristic  number  with  many  elementary  divisors,  and  m  has  many  non-zero 
components,  this  leads  to  cumbersome  calculations..  If  we  replace  m  by  in, 
with  only  one  non-zero  component  for  each  distinct  characteristic  number  X ,  the 
work  is  easier.  As  an  example,  suppose  that  D  is 
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If  m  =  0 ,  m 
1  2 


m^  =  0,  then  (14)  becomes 


-1 


-1 


h>6x  ( 6 1/4X  +Cl/2)  +  o  2  ( 


-&/2X  -  5  x/ 4X 


3  +  c2/2)2  +  63“1(c3/2)2 


Clearly  c3  is  unrestricted.  If  both  c  and  c2  are  non-positive,  the  right  side 

can  be  made  arbitrarily  small  for  suitable  choice  of  A.  if  c^>  0,  then  for  any 

2  — 

5^,  the  second  term  is  minimal  if  b^  =  b  2k  -  kc^  >  0  with  minimal  value 

4  _ 

6  ^/2X  -  c^/k  .  Thus  sufficient  for  absolute  stability  is 


h  >  96^/l6\4  +  c12/461  +  c^/4kZ  -  c^/k  . 


2  _ 

If  6^  =  2X  |  c  l/3  the  right  side  is  again  minimized  and  a  stability  criterion  can 
be  written  as 

< 


h  >  €  (c2)  { 


3|C1I+C1  °2 


. 


4X 


where  s  ( t)  =  0  if  t  <  0 ,  while  €  =  1  if  t  >  0 . 

Suppose  now  that  the  reduction  of  ( 1)  to  (  2)  had  yielded  =  m2  =  m3  =  1  • 


Then  U  can  be  taken  as 
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1 

0 

0 


1  ' 

1 

1 


0 

0 

0 


since  c  =  U  c,  we  have 


c  =  c 
1  1 


C2  =  °1  +  °2  +  C3 


and  sufficient  for  absolute  stability  is 


h  >  €  (c1  +  c2  +  c3) 


In  the  present  case  the  saving  in  effort  is  small  but  noticeable.  If  the  exponents 
of  the  various  elementary  divisors  had  been  higher,  the  saving  would  be  more 


significant, 
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APPENDIX 


Kriuchkov  [9]  (in  our  notation)  represents  the  fundamental  matrix 
At 

e  of  the  homogeneous  system  yl  =  Ay  as  a  linear  combination  of  the  terms  of 

•  \  t 
]  k 

the  form  t  e  with  matrix  coefficients  which  are  the  so-called  components  of 
A.  In  this  form,  the  highest  power  t"1  which  appears  is  one  less  than  the  power 
of  (X -X  )  which  appears  in  the  minimal  polynomial  of  A(cf.  e.  g. ,  [5],  p.  110  ff. ) . 
If  k  is  the  degree  of  the  minimal  polynomial,  then  the  solution  of  the  system 


(16) 


'  y'  =  Ay  +  up 

< 

u  —  f( s)  ,  s  =  (b,  y) 


is  represented,  using  the  variation  of  constants  formula,  as 

At  .  «  £ 
y  =  e  y  +  G§ 
o 


where  G  is  a  certain  constant  n  x  k  matrix  and  £  is  a  particular  solution  of 


(17) 


=  A  £  +  mf(  s)  . 


Here  A  is  a  k  x  k  matrix  in  Jordan  normal  form  whose  characteristic  ( and 
minimal)  polynomial  is  the  minimum  polynomial  of  A,  while  the  projection  of 
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m  on  each  invariant  subspace  of  A  has  component  one  along  the  characteristic 
vector  and  zero  elsewhere.  Should  X.  be  a  characteristic  root  with  positive 
real  part,  of  greater  multiplicity  in  the  cnaracteristic  equation  than  in  the 
minimal  equation  of  A  the  stability  of  (.17)  does  not  imply  the  stability  of  ( 16). 
This,  however,  is  not  the  Lur’e  transformation,  as  constructed  above.  The 
calculations  above  show  that  there  is  a  vector  m,  and  nonsingular  matrices 
B  and  V  such  that  y  s  VBz  transforms  ( 16)  into 

'  z1  ~  Dz  +  mf(  s) 

< 

s  —  { c ,  z) 

-1  T  T 

where  D=B  AB,  c  =  B  V  b,  and  V  is  a  matrix  which  commutes  with  A. 

In  this  instance,  equation  { 5)  becomes 

( 5” )  VBm  »  p. 

Setting 


V  =  BUB 


3 


we  obtain 
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(6') 


Um  =  p 


where  U  is  an  arbitrary  matrix  commuting  with  D .  It  is  easy  to  see  that 
the  construction  used  for  DU  is  now  a  construction  for  U.  This  Lur' e 
transformation,  however,  clearly  preserves  stability  and  instability  properties. 

Acknowledgment:  I  wish  to  thank  Dr.  P.  C.  Parks  for  his  valuable  comments 


on  an  earlier  version  of  this  paper. 


-30- 


#392 


BIBLIOGRAPHY 

[  lj  Lur’e,  A.  I.  ,  Some  Nonlinear  Problems  in  the  Theory  of  Automatic  Control. 
Gos.  Isdat.  Tekh.  Theor.  Lit.  Moscow  ( 195 1)  Available  in  English  from 
Her  Majesty’s  Stationery  Office,  London  (1957). 

[  2]  Letov,  A.  M.  ,  Stability  in  Nonlinear  Control  Systems.  Princeton  University 
Press,  Princeton,  New  Jersey,  pp.  33-47. 

[3]  La  Salle,  J.  and  Lefschetz,  S.  ,  Stability  by  Liapunov’s  Direct  Method. 
Academic  Press,  New  York  (1961). 

[4]  Gibson,  J.  E. ,  McVey,  E.  S. ,  Leedham,  C.  D. ,  Rekasius,  Z,  V., 

Schultz,  D.  G. ,  and  Sridher,  R. ,  Stability  of  Nonlinear  Control  Systems  by 
the  second  method  of  Liapunov.  AFMDC -TR-61-6,  Vol.  II,  U.  S.  Dept,  of 
Commerce,  Office  of  Technical  Services,  Washington  25,  D.  C.  (1961). 

[5]  Hahn,  W. ,  Theorie  und  Anwendung  der  Direkten  Methode  von  Liapunov. 
Springer,  Berlin  (  1959) . 

[6]  Turnbull,  H.  W. ,  and  Aitken,  A.  C.,  The  Theory  of  Canonical  Matrices.  ... 
Blackie  and  Sons,  London  (Second  Edition),  (1948) . 

[7]  Gantmacher,  F.  R. ,  The  Theory  of  Matrices,  Vol.  I.  Chelsea  Publishing  Co. , 
New  York  (  1959) . 

[8]  Krasovskii,  N.  N. ,  Stability  of  Motion.  Stanford  University  Press,  Stanford, 
California  (  1963) . 

[  9]  Kriuchkov,  N.  B. ,  Investigation  of  a  Pecularity  of  Nonlinear  Equations  of 
Control  Systems  in  the  Case  of  Multiple  Roots  by  Application  of  the  Theory 
of  Matrices,  PMM  Vol.  24,  No.  4,  I960,  pp.  749 -754(  translated  in  Applied 
Mathematics  and  Mech.  Vol.  24,  No.  4,  1960  pp.  1125-1134.) 


